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Swimming cells often have to self-propel through fluids displaying non-Newtonian rheology. While 
past theoretical work seems to indicate that stresses arising from complex fluids should systemat¬ 
ically hinder low-Reynolds number locomotion, experimental observations suggest that locomotion 
enhancement is possible. In this paper we propose a physical mechanism for locomotion enhance¬ 
ment of microscopic swimmers in a complex fluid. It is based on the fact that micro-structured fluids 
will generically phase-separate near surfaces, leading to the presence of low-viscosity layers which 
promote slip and decrease viscous friction near the surface of the swimmer. We use two models to 
address the consequence of this phase separation: a nonzero apparent slip length for the fluid and 
then an explicit modeling of the change of viscosity in a thin layer near the swimmer. Consider¬ 
ing two canonical setups for low-Reynolds number locomotion, namely the waving locomotion of 
a two-dimensional sheet and that of a three-dimensional filament, we show that phase-separation 
systematically increases the locomotion speeds, possibly by orders of magnitude. We close by con¬ 
fronting our predictions with recent experimental results. 


I. INTRODUCTION 

Over the past few decades, problems on life at low Reynolds number have received significant attention, from both 
the physics and biological communities [1-5]. For the most part, theoretical and experimental studies have focused 
on cell locomotion in Newtonian fluids, with an emphasis on the interplay between biological actuation and whole¬ 
cell response, for example the relationship between cell geometry, waving actuation, and the resulting swimming 
kinematics [4, 6]. 

Studies on cell motility in fluids typically focus on one of four types of cells bacteria [7], spermatozoa [8], ciliates 
[9], and single-celled planktonic cells [10] - while recent effort considered larger, multi-cellular organisms [11]. In all 
cases, the biological environments that the microorganism encounter may be Theologically complex. For example, 
Helicobacter pylori , a bacterium that causes inflammation in the stomach, swims through gastric mucus to protect 
itself from the acidic environment [12]. Mammalian spermatozoa have to progress through highly-elastic cervical 
mucus, an important phase in reproductive process [13]. 

Extending our understanding of cell locomotion in Newtonian fluids to complex, gel-like or viscoelastic environments 
is a nontrivial task. One quintessential question, whether non-Newtonian stresses in a complex fluid help increase or 
decrease the swimming speed of the cell, remains in many ways an open problem. 

For example, bacteria such as Leptospira and Escherichia coli , swim more rapidly in gel-like unbranched polymer 
solutions than in Newtonian fluids [14]. In contrast, for the nematode Caenorhabditis elegans undergoing undulatory 
swimming, the speed was observed to decrease in a slightly shear-thinning polymeric fluid with strong elastic stresses 
[15]. Similar disparities are observed experimentally for bio-inspired synthetic swimmers. Force-free rotating helices 
show a transition from hindered to enhanced swimming in constant-viscosity Boger fluids [16] while a cylindrical 
version of Taylor’s swimming sheet displays both increase and decrease as a function of the rheology of the fluid [17]. 
In contrast, externally-actuated flexible-tail swimmers show a systematic increase of locomotion speeds in viscoelastic 
fluids [18]. 
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Various numerical and theoretical studies have also addressed this problem, focusing on viscoelastic fluids follow¬ 
ing Oldroyd-B rheology. Small-amplitude asymptotic studies for waving swimmers with fixed shapes predicted a 
systematic decrease of swimming velocity [19-21]. Subsequent numerical work for finite waving sheets with large 
tail amplitude showed that an increase was possible for order one Deborah numbers [22]. Numerical simulations 
following the helical experiments in Ref. [16] confirmed the transition from slow small-amplitude swimming to fast 
large-amplitude locomotion [23]. Integral theorems for small-amplitude motion showed that the superposition of 
multiple waves could also lead to a enhancement transition for a range of Deborah numbers [24]. 

We thus see that theoretical, computational, and experimental studies showing both increases and decreases have 
been put forward, and the challenge is now to rigorously untangle the various physical (and sometimes, biological) 
effects. In particular, while we now understand how viscoelastic stresses are able to decrease swimming speeds, 
physical mechanisms leading to locomotion enhancement are less clear. Recently, the flexibility of the swimmer in 
response to complex stresses was shown to allow for an increase in the swimming speed [25, 26]. In this paper, we 
propose a different physical origin for the observed swimming enhancement. Instead of focusing on the new non- 
Newtonian stresses in the fluid, we address one of the consequences of having a structured fluid, namely the fact that 
it is expected to phase-separate near the body of the swimmer, leading to the well-known phenomenon of apparent 
slip. 

For a variety of complex fluids with a microstructure dispersed in a solvent, in particular polymeric fluids and 
suspensions, the presence of a boundary leads to static phase separation at equilibrium: the concentration of the 
solute drops near the wall which is covered instead by a thin solvent layer [27]. In the case of rigid suspensions, 
purely excluded-volume interactions lead to solvent-rich regions near the surface, and the effect is larger for Brownian 
particles for which the presence of a wall breaks the geometrical isotropy [27]. In the case of polymers, random coils 
would be distorted if too close to the wall, and thus they are driven by entropy away from the boundary. 

In all cases, the solvent-rich fluid near the surface has a viscosity much smaller than that of the bulk fluid. As seen 
in many situations, in particular flows in capillary tubes and in porous media [28], this difference in viscosity leads to 
apparent slip when a flow is set up, which is best illustrated in the case of a shear flow (Fig. lb): If shear is imposed 
in a fluid with a thin-viscosity layer, the difference in viscosities will lead to a difference in shear rates, and as a result 
the flow in the high-viscosity bulk will not extrapolate to zero on the solid surface, but below it, indicating an overall 
decrease of stresses acting on the surface. The fictitious distance below the surface where the fluid velocity in the top 
fluid goes to zero is the (positive) apparent slip length. Microscopically, the no-slip condition is of course not violated, 
but given that the typical thickness of the solvent layer is much smaller than the other, macroscopic length scales in 
the problem of interest, the no-slip boundary condition appears not to hold for the bulk fluid. 

There are two classical ways to theoretically model apparent slip in complex fluids. The first model is to sim¬ 
ply replace the no-slip boundary condition on the surface by one which the tangential velocity is allowed to slip. 
Experimentally-measured slip length has been shown to depend, sometimes in a complex manner, on the shear stress 
at the wall [28-31]. The assumption usually done is to adopt Navier’s slip length model [32] and assume that the 
slip velocity at the wall is linearly proportional to the wall shear rate, with a proportionality constant with unit of 
length, called the slip length, and which we will denote A in this work [33]. As noted above, the slip length measures 
the (fictitious) distance below the boundary where the velocity would extrapolate to zero, and it is zero in the case 
of a no-slip boundary (Fig. la). A second procedure to model phase separation is to explicitly assume the presence 
of two fluid layers. The top layer, semi-infinite, has bulk viscosity /zi, while the bottom layer near the surface has a 
finite thickness h and lower viscosity H 2 < Hi (Fig. lb). For a shear flow in this unidirectional setup, the velocity in 
the thin layer satisfies the no-slip boundary condition while that the flow in the bulk fluid extrapolates to zero at the 
equivalent apparent slip length A = (hi/h? — 1 )h. 

In this paper we consider these two different physical models of phase separation and investigate their consequences 
on waving locomotion. In §11 we first examine the model with a finite apparent slip length, and apply it to two 
canonical setups for low-Reynolds number locomotion, namely the small-amplitude swimming of a flexible sheet [1] 
and that of a flexible filament [34]. In §111 we then consider the same two setups in the situation where the phase- 
separated fluid is modeled as a two-fluid layer. In all cases we are able to derive the swimming speed for each swimmer 
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(a) 




FIG. 1: (Color online) The two models of apparent slip due to phase separation considered in this paper; (a): a single-phase 
continuum fluid with a finite apparent slip length A; (b): a two-fluid domain with viscosity pi in the bulk and a thin low-viscosity 
layer of thickness h and viscosity p 2 < pi- 


analytically and we compare it to the case for the homogenous Newtonian fluid with a no-slip boundary condition. 
We demonstrate that the phase separation leads to a systematic enhancement of the swimming speed, and suggest 
that this might play a role in the recently-measured swimming enhancement at low-Reynolds numbers. 


II. SWIMMING IN A FLUID WITH FINITE APPARENT SLIP LENGTH 

In the first section we assume that the phase separation in the fluid can be adequately captured by an effective slip 
length A acting on a Newtonian fluid satisfying the Stokes equations 

'S/p = /rV 2 u, V ■ u = 0. (1) 

On a fluid-solid boundary S, the jump in normal velocity is zero by mass conservation while the jump in tangential 
velocity is proportional to the local shear rate. If the velocity in the fluid is denoted u, these boundary conditions 
can be mathematically expressed as 


[ n ' u ] I s = °> ( 2a ) 

[n x u] | s = 2A(n x (E ■ n))| s , (2b) 

where [...] is used to denote a jump, n is the normal to the boundary, E is the symmetric rate-of-strain tensor (i.e. the 
symmetric part of the velocity gradient tensor), and A the slip length. 


A. Two dimensional waving sheet 

We first consider a two-dimensional swimmer in the form of flexible sheet self-propelling in the fluid by passing 
waves of normal deformation. This is the classical setup originally proposed by Taylor [1], and the material points of 
the sheet, (x s ,y s ), are assumed to vary in space and time as a simple traveling wave of deformation 

y s = bsink(x — ct), x s = x, (3) 

where b is the wave amplitude, k the wave number and c the wave speed along the x direction (see notation in 
Fig. 2). We solve the problem assuming that the amplitude is small compared to the wavelength, and thus consider 
the asymptotic limit where e = bk is a small dimensionless number. Nondimentionalizing the equations using fc -1 as 
relevant length and w -1 = (fcc)^ 1 as intrinsic time scale, the wave deformation becomes 


y s = esin(x — t) = esin£, 


(4) 
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FIG. 2: (Color online) Geometry of the two-dimensional waving sheet setup in a fluid with finite slip length. Here b denotes 
the small waving amplitude, A the wavelength, and c the wave speed along the x direction. The sheet is assumed to swim with 
speed U in the negative x direction. The small circle zooms in on a portion of the swimmer surface and illustrates the presence 
of a finite slip length A. 


and for convenience we drop the “bars” in what follows. 

The velocity of material points on the sheet is thus given by 

u s = (0, —e cos £). 

Applying Eq. (2), on the surface of the swimmer we have 


( n ' u ) 


Os,y») 


= n ■ u. 


(“ X U ) | = 2A ( nX ( E ' n ))|(x a ,y a ) +nXU *’ 


( 5 ) 

(6a) 

(6b) 


with a normal vector n explicitly given as (1 + e 2 cos 2 £) _ s [—ecos£, 1], and where A = kA is the nondimensionalized 
slip length. Here we keep the “bar” notation for the slip length to ensure that the final result is formally dimensionless. 

In order to solve for the Stokes equations in the fluid, Eq. (1), we employ a streamfunction ip such that the velocity 

d'ljj d'tjj 

components of u = \u. u] are defined as u = —— and v = — -7—. The boundary condition, Eq. (6), then becomes 

ay ox 


.dip . dip, 

ecos£ —- (x, esm£) + —— (x, esm( = ecos£, 
ay ox 


.dip 


dip 


ecos£ — (x, esin£)—— (:r,esin£) = A(l + e cos £) 2 < (1 — e“ cos £) 


1 dx 


dy 


(7a) 

(7b) 


<9 2 Vg, . ^ . , 

alaU.esmO - ^(£,esm£) 


dx 2 


d 2 ip 


+4e cos p Q x Q y (Ci e sin £) + e 2 cos 2 £. 

The Stokes equation, Eq. (1), transforms into the biharmonic equation for ip [35] 

v 4 v> = o. 

In order to obtain the asymptotic solution for the swimming velocity, we expand the streamfunction in powers of e by 

ip = eip^ + e 2 ip^ 2) + e 3 ip( 3) + ... (9) 


( 8 ) 


Denoting the velocity of the swimming sheet as —Ue x in a quiescent fluid, we move in the swimming frame and thus 
the velocity at infinity is given by u (y —> oo) = Ue x . 

Since ip satisfies the biharmonic equation, and is equal to Uy at infinity, we construct the general solution as [36] 


ip w = V^ + UWy, (10a) 

V> (2) = v} 2) + E 2 (2) + t/ (2) y, (10b) 
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where 


V n = (A n + B n y)e~ ny sin n£ + {C n + D n y)e~ ny cos n£. 
At first order in e, Eq. (7) becomes 

dipW 


dx 

dV> (1) 

dy 


(a;,0) = cos 

(x,0) = A t82 * m d2 *" 


\ dx 2 


dy 2 

J\ 

Or,0) 

b[ 1] 

= 1, C\ 


(ii) 

(12a) 

(12b) 


streamfunction at first order is 

ipA) = (1 + y)e~ v sin£, (13) 

which is the same as the no-slip case. This can be rationalized by notating that the first-order shear rate is given by 


/ 9 2 V> (1) < 9 2 V> (1) \ 

\ dx 2 dy 2 ) 


(*,o) 


= (—2 ye y sin^) 


= 0, 


(14) 


O,o) 


which makes the problem equivalent to the no-slip case. 

As no propulsion occurs at order e, one needs to carry the calculation to order two in order to obtain the leading- 
order swimming speed. At order e 2 , the boundary conditions are 


Qxpd) dipt 2 ' 1 d 2 ip W \ 

cos £ H7 I HZ ^ sin £' 


dy 


dx 


dxdy J 


dipdipt 2 ') d 2 ip W \ 

cos£—- - -sin£- 


dx 


dy 


dy 2 J 


O,o) 


O,o) 


0, 

A 


Ud 2 ipW 

d 2 ip( 2 ^ \ 

\ \ dx 2 

dy 2 ) 


+ 


(15a) 

(15b) 


( d 3 ip^ d 3, ip^ \ 

Sm? \dx 2 dy~^dy^) 


+4 cos £ 


d 2 ip^ \ 


O,o) 


1 1 
2 + 2 


cos 2£. 


dxdy J 

Substituting the expansions for the streamfunction into this condition, we obtain 

U {2) = ^ + A. 

Comparing this result with the no-slip case, and coming back to the dimensional variables we finally have 

[/(2) 


u: 


( 2 ) 


= 1 + 2kA. 


(16) 


(17) 


no—slip 


Since the slip length is always positive, we obtain in this first situation that the swimming speed is always enhanced 
by apparent slip. 


B. Three-dimensional waving filament 

We now apply the same apparent-slip model to the case of a three-dimensional waving filament, the geometry of 
which is shown in Fig. 3 [34]. We consider a cylindrical filament of radius p deforming as a traveling wave in the (x, z) 
plane where 2 is along the filament axis and x is perpendicular to it. We denote by <5 the amplitude of the filament 
deformation in the x direction. 
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FIG. 3: (Color online) Deformation with amplitude 5 in the circular cross-section of a three-dimensional waving filament of 
radius p. The blue solid line represents the current position of the filament and the gray dashed line the average location of 
the cross section. 

The location of the material points on the surface of the filament, using cartesian coordinate centered on the average 
location of the cylinder (see Fig. 3), is given by 

r s = (d + pcosG)x + psind y + z z, 6 = bsmk(z + ct), (18) 

where <5 is the distribution of waving amplitudes along z direction, p the filament radius, k the wave number, and c 
the wave speed. Here again we assume that the amplitude is small compared to the wavelength, and compute the 
result in the limit where e = bk is small. 

Nondimentionalizing the equations similarly to the previous section, we have the surface of cylindrical filament 
described as 


r s = (5 + pcos 6 »)x + psin 0 y + z z, S = esin (,3 + t) = esins, 


(19) 


and here again we keep the “bar” notation for the dimensionless radius, p = kp. Similarly to the two-dimensional 
case, the velocity on the boundary satisfies the conditions 


(n-u)|,,,_ ecosscosd, ( 20 a) 

v ' I (d+pcos0,psin0,2:) 5 v y 

(n x u)L,, _ = 2A{n x (E • n)}|, x . _ esinflcoss z. (20b) 

v ' I (d+pcosQ,psinQ,z) L v I (d+pcos 0,psin 6,z) v J 

Since the inextensibility condition contributes to the dynamics at orders higher than two [34], the vector normal 
to the surface is e r at the order relevant for this calculation. We can then expand the velocity around the average 
position as 


u(<5 + p cos d, p sin 0, z) 


[u + x ■ (Vu)esins] 
du 


u + e sin s 


dr 


r—p 

cos 8 — 


/1 du 
\rd0 



+ 


v + e sin s 


dv 

dr 


cos v — 



G r 

r—p 

ee 

r=p 


+ 


w + e sin s 


dw 1 dw . 

— cos - 7777 sm 9 

dr r dO 


e*. 

r=p 


Expanding the velocity in the fluid, u, asymptotically in powers of e, 


( 21 ) 


u=eu (D +e 2 u (2) + .. 


(22) 
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and substituting the expansion and Eq. (21) into the boundary condition, we obtain that at first order, 


■ W (p,0,z) 

= cos s cos 0. 


m (p,0,z) 

— fdv (1) 

ul 1 ) 

(dr 

r 

m (p,0,z) 

— ( dw^ 

du-W ' 

= A 

+ 


(r=p) 


— sin U cos s, 


(r—p) 


Solving for the fluid velocity using separation of variables we get 

= u q {r) cos0coss, = v q (r) sin0coss, = w q (r) cos 9 sins, 

while the first-order boundary conditions are simplified to 


Uq(P ) = 1, 

Vq{P ) = A 


V'q(p) - 


V q (p) 1 


- 1 . 


j(p) = A K (p) - x ] • 


(23a) 

(23b) 

(23c) 

(24) 

(25a) 

(25b) 

(25c) 


Similarly to the original problem treated by Taylor in the case of a no-slip filament [34], the radial dependence of 
the velocity is given by a combination of modified Bessel functions as 


u q (r) = BI\ 2 (r) + CK 0 (r) + ArKi(r)+EI 2 {r) + FI 0 (r) + Drli(r), 
v q (r) = BK 2 (r) — CK 0 (r)+EI 2 (r) — FI 0 {r), 
w q (r) = BK x (r) + CK 1 (r) + A[rK 0 (r) - K 1 (r)]-EI 1 (r) - FIi(r) 
—D[rIo(r) - h(r)], 


(26a) 

(26b) 

(26c) 


with a minus sign in Eq. 26c coming from different properties between solution of the first and second kind. For 
the boundary conditions at infinity, we have D = E = F = 0 and the other three unknown constants A 1 B, C can 
be obtained by plugging Eq. (26) into Eq. (25). Writing C as C = C nu /Cd e , we obtain the following lengthy (but 
analytical) expressions 


C™ = pK^p) 2 - 2pK 0 (p)K 2 (p) + 2K 1 (p)K 2 (p) + A [-^Ko^K^p) 
—K 0 (p)K 2 {p) + 2K 1 (p) 2 + - -pj K 1 (p)K 2 (p) + K 2 {pY 


(27a) 


+A 2 


2K 0 (p)K 1 (p) - 2pKYp) 2 + -_K 0 (p)K 2 (p) 


-6 Kd{p)K 2 {p) + -_K 2 {p) 2 
P 

C de = -2pK 0 (p) 2 K 2 (p) + 2K 0 (p)K 1 (p)K 2 (p) + pI^K^p) 2 + pK^pfK^p) 


(27b) 


+A 


- , -K 0 (p) 2 K 1 (p + 2pK 1 (p) 3 + 2K 0 (p)K 1 (p) 2 + K 0 (p)K 2 (p) 2 


-K 0 (p) 2 K 2 (p) +4K 1 (p) 2 K 2 (p) + Z Kl (p)K 2 (p) 2 


-2p)K 0 (p)K 1 (p)K 2 (p) 
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A 2 ^KYpfKYP) + -K 0 {p) 2 K 2 {p) 
P 


+ -Ko(p)K 2 (p) 2 + 2R\(p)K 2 (p) 2 
P 


B = -^{-pKiip? +A 


~^pK 0 (p)K 1 {p) - ^Ki(p)K 2 (p) 


— 2K 2 pKi(py 


(28) 



and 


A=—{2K 1 (p)K 2 (p) + A 

*-'de 


2K\(p) 2 + 3K 0 (p)K 2 (p) + -A h{p)K 2 {p) + K 2 (p) 2 

P 


+A 2 1 2Ko(p)Ki(p) + -_K 0 (p)K 2 (p) + 2K\{p)K 2 (fi) + -I< 2 (p) 2 
P P 


(29) 


According to Eq. (24), the time-averaged swimming speed is zero, and as expected we thus need to consider the 
problem at order e 2 . 

Following the slip boundary conditions, the ^-component of the flow at order two, w^ 2 \ satisfies on the boundary 


+ w' cos 2 9 sin 2 s H— w a sin 2 9 sin 2 s 


=A 


r—p 


dw^ du^ 
dr dz 


+ sin 2 s\(w” — u' ) cos 2 9 


lCy~ w, q + u 4 + v i) sin2 6 


(30) 


o I / 2/i ^q + Vq , n 

+ cos 2 s | u' cos 2 9 H— -- sm ( 

H r 


r—p 


Averaging this equation in time and along the azimuthal direction, we obtain explicitly the swimming velocity, U^ 2 \ 
as 

\ (w„ 


4 [A 2 ) = w'Jp) + - A 

P 


w q - - 
. P \ P 


-Ki(p) - K 2(p) - Ai^i(p) 
P 


B 


-K 1 (p) - K 2 {p ) - AK 1 (p) 
P 


C 


(31) 


+ <{ 2K 0 (p) - pKi(p) - -K^p) + K 2 (p) + A[5Ad(p) - pK 2 (p)\ } A. 


The result in Eq. (31) can also be evaluated in the no-slip case by simply setting A = 0, and we recover Taylor’s 
result, namely 


U, 


( 2 ) 


pKi(p) 2 K 2 {p) - pK 0 (p)K 2 (pf 


" sli P -2pA 0 (p)Ad(p) 2 + 4pK 0 (p) 2 K 2 (p) - 4A” 0 (p)AA(p)AT 2 (p) - 2pK 1 {pYK 2 {p) 


(32) 


The ratio between the swimming speed in the slip case to that in the no-slip situation, /U^_ slip , is plotted in 
Fig. 4 as a function of the dimensionless filament radius (Fig. 4a) and the dimensionless wave number (Fig. 4b). As in 
the two-dimensional situation, the presence of slip is seen to always lead to faster swimming than in the no-slip case, 
and here the effect can be potentially very large (the applicability of these results to recent experiments is discussed in 
§IV). In Fig. 4 we see that the swimming speed increases monotonically when either the length scale of the swimmer 
cross section ( p ) or the typical length scale of the waving motion, k~ 1 becomes smaller than the slip length. In the 
opposite limit, the no-slip result, Eq. (32), is recovered when all length scales are much larger than A. 

We further note that it is possible to compute the swimming speed for small slip length as a power expansion in 
kA (i.e. A), CA 2 ) = £/q 2) + ( kA)uf’' > + ..., with being the no-slip swimming speed. When kp (i.e. p) increases to 
infinity, the radius of the cylinder becomes much larger than any other length scale, and we recover /Uq 2 ' 1 = 2, 
leading to U ( - 2 ' 1 /U^_ slip = 1 + 2/cA, which as expected agrees with the results for the two-dimensional sheet. 

We conclude by pointing out that although the calculation above was carried out in the case of planar waving 
deformation, similar algebra would govern swimming by propagating helical waves [37], and in that case p would 
be the radius of the helical flagellum, or that of the bundle of flagella in the case of bacteria with multiple flagellar 
filaments such as E. coli. 
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FIG. 4: (Color online) Ratio of the swimming velocity in the case of slip to the no-slip value, /U^_ Alp . (a): Dependence on 
the dimensionless filament radius, kp, for three wave numbers (kA = 0.01, 0.1 and 0.5); (b): Dependence on the dimensionless 
wave number (with kp = 0.01, 0.1 and 0.5). 


III. SWIMMING IN A TWO-FLUID DOMAIN 

In the previous section, we modeled the influence of phase separation as due to a finite apparent slip length, 
and showed that it leads to a systematic enhancement of the swimming speeds. In order to provide an alternative 
microscopic physical picture, we instead consider in this section a second model where we include explicitly the 
presence of a low-viscosity layer near the surface of the swimmer, and we apply it to the two canonical swimmers 
(waving sheet and filament) considered in the previous section. 


A. Two dimensional waving sheet 

We first consider the swimming sheet setup shown in Fig. 5. The fluid is composed of two domains: the bulk fluid 
has viscosity p\ while the thin layer near the swimmer has a smaller viscosity p 2 (pi > ^ 2 )- The average distance 
between the sheet and the fluid-fluid interface, which is the thickness of the low-viscosity layer, is denoted h and 
assumed to remain constant (this assumption is discussed in §IV). All other notation are similar to the ones in §11 A. 


Ml, ^1 



FIG. 5: (Color online) Geometry of two-dimensional waving sheet swimming in a two-fluid domain. The bulk fluid has viscosity 
pi and the thin layer near the swimmer, of mean thickness h, has viscosity p 2 < Pi- The two streamfunctions in the fluids are 
denoted ip i and i/’ 2 - The interface between the two fluids is assumed to remain flat. 
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Following the same nondimentionalization, we now have the no-slip boundary conditions on the sheet 

— (x,esm£) = 0, — (x,esm£) = -ecos£, 
while at inhnity we have the unknown swimming speed 


oipi, s TT clip i 

oo) = —U, oo) = 0. (34) 

At the flat interface between the two fluids, we have continuity of velocities 

w (x ' S) = W 1 *’ S)> lt {x ’ 5) = S) = (35) 

and here again we keep the “bar” notation for the thickness h = kh. Together with continuity of tangential stresses, 
which is written as 


a 2 ip l a 2 v>i 


d 2 ^2 9 2 *l>2 


\ dy 2 dx 2 ) \ {x - h) V 9y 2 9x 2 ) 

where (3 denotes the ratio of viscosities, /3 = H 2 / Hi < 1- 

As in §11A we solve the problem as a perturbation expansion in e. The general periodic solution for to the biharmonic 
equation which vanishes at inhnity is obtained by separation of variables as 


which we use to expand 1 pi as 


V n = (a n y+ b n )smn£e ny + (c n y + d n ) cosn£e ny , 


4>i = + e 2 V>2 2) + • • • 


ip[ m) = -U {m) y + Vf m) + ■ ■ ■ + 


Similarly, in the second domain we expand if )2 as 

V>2 = eV’2 1 ’* + + ■ ■ • (40) 

^rn) = w (jn) + . . . + W (rn) + ^m ) y 2 + (41) 

where W n is the general periodic solution for to the biharmonic equation in a finite domain obtained by separation 
of variables 

W n = [( A n y + B n ) sin ni ; + ( C n y + D n ) cos n£] sinh ny 

( J 

+ \{E n y + F n ) sin + (G n y + H n ) cos n£] cosh ny. 


Expanding the boundary conditions around y = 0, we obtain at first order 


(on sheet) 


^-(i,0) =0, ?h-fa0) = -cos( 


(on interface) —— (x, h ) = —— (x, h ), —— (x, h ) = —— (x,h) = 0 
oy oy ox ox 


d 2 ip{ 1) d 2 ip{ 1) 

dy 2 dx 2 


d 2 %d 2 il) ^ 
dy 2 dx 2 


(at inhnity) 


~^—(x,oo) = -C/ (1) , -^-(*, 00 ) = 0. 
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FIG. 6: (Color online) Ratio between the swimming velocity in the two-fluid domain, U^ 2 \ and the velocity in the presence of 
a single fluid, Us, o, as a function of the dimensionless distance between the swimmer and the interface, kh , for three values of 
the viscosity ratio: /? = 0, 0.1 and 0.5. 


Substituting these boundary conditions into the general solution we obtain the coefficients at first order 

t/ (1) = X W = j/ 1 ) = C W = rfW = c[ 1] = d[ 1) = 0, 


,(i) 


4 1} = 


A™ = 


b[ 1} = 


/3he h sinh h 


sinh 2 h — h 2 + /?(sinh h cosh h — h) 


—(3h 2 e h sinh h 


sinh“ h — h 2 + /3(sinh h cosh h — h) 
sinh 2 h + /3(sinh h cosh h) 
sinh 2 h — h 2 + /3(sinh h cosh h — h)' 
sinh h cosh h + h + f3 cosh 2 h 
sinh 2 h — h 2 + /3(sinh h cosh h—h) 


At second order, the boundary conditions become 


(on sheet) 
(on interface) 


dy 


(x, 0 ) + sin£ 


d 2 ^ 1] 

dy 2 


(z,0) = 0, 


dy 

d 2 ^^ d 2 i/j^ ^ 

dy 2 dx 2 , 


d ^\ x - h ) = d $^{xrh). 


dy 


= P 


' 5 2 4 2) d 2 ^ 2) \ 


(x,h) 


dy 


dx 2 


( x,h ) 


leading to the second-order swimming speed as 


J7< 2 > = = Af) _ 


and therefore 


uw = I + — 2 _ _ h2 + P h _ - 

2 sinh h — h 2 + /?(sinh h cosh h—h) 


(44a) 

(44b) 

(44c) 

(44d) 

(44e) 

(45a) 

(45b) 

(45c) 

(46) 

(47) 


( 2 ) — 

In Eq. (47), the first term is the one-fluid classical result of Taylor (Uso = 1/2, recovered when h —> oo) and the 
second fraction is always positive since h > 0. As a consequence the swimming speed for a waving sheet in a two-fluid 
domain is always faster than in a homogeneous Newtonian fluid. 
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FIG. 7: (Color online) Waving motion of amplitude S in the circular cross-section of a three-dimensional waving filament of 
radius p in a two-fluid domain. The solid blue line (inside) indicates the instantaneous position of the filament, the dashed gray 
line the average location of the cross section, while the solid red line (outside) shows the interface between the thin low-viscosity 
layer (mean thickness h; dynamic viscosity p 2 ) and the bulk fluid (viscosity pi > pi). 

We display in Fig. 6 the ratio /U& as a function of the dimensionless distance to the interface, kh (i.e. h). 
We observe that the increase of the swimming speed can become very large when the thickness of the low-shear layer 
is smaller than the wavelength of the swimmer. We also see that the overall conclusions and speed ratios are rather 
insensitive to the exact value of the viscosity ratio, /3. 


B. Three dimensional waving filament 


In this final section, we extend the two-fluid scenario to the case of three-dimensional waving filaments. In this 
case, the geometry of the cross-section, shown in Fig. 7, is analogous to the one addressed in §IIB with the added 
ingredient that we now have two fluids. The thin, low-viscosity layer, has mean thickness h and dynamic viscosity P 2 
while the bulk has viscosity pi > p 2 - All other notation are similar to the ones used in §IIB. 

On the surface of the filament (dimensionless form), r s = (6 + pcosd)x + p sin fly + zz, we have the distribution of 
surface velocities U2(<5 + pcos0,9,z) = ecossx, which can be written in polar coordinates as 

112(8 + pcosO, 9, z) = ecos^coss, v 2 (S + pcos0,9,z) = -esintfcoss. (48) 


Expanding the velocity components around the averaged position of the surface Eq. (21), we obtain at first order 

(p, 0, z) = cos 9 cos s, v^\p, 0,z) = — sin 0 cos s, w^\p, 9, z) = 0. (49) 

We assume that the interface r = h + p = l undergoes no radial motion and apply continuity of the tangential 
components of velocities and traction leading to the conditions 


u\ J (Z,6»,z) = u\ (l, 9, z) = 0, 
v[ 1] (i,o,z) = vi 1] (i,e,z), w { i\i,e,z) = wi 1] (i,e,z), 


dw\ 


( 1 ) 


du 


(i)' 


dr 


dz 


” = l 




(i) 


du 


(i)' 


dz 


d(v^/r) 1 chij 1 ) 

= P 

d(^2 ) / r ) 

dr r dd 

r—l 

dr r dd 


" = l 


r—l 


(50a) 

(50b) 

(50c) 

(50d) 


where, as in §111 A, /3 denotes the ratio of viscosity, /? = p 2 /Pi < 1- 
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Using separation of variables, we write 

= u\ q (r) cosllcoss, = ui 9 (r) sin 0 cos s, = wi q (r) cos 9 sins, (51a) 

= u 2q (r) cos 9 coss, = v 2q (r) sin0coss, = w 2q {r) cos 9 sins, (51b) 

and substituting Eq. (51) into the Eq. (50), the boundary conditions become 

U2q(p) = 1 , V 2q (p)=-1, w 2q {p) = 0 , 

Ul q (l) = U 2 q{l) = 0 , V\q(l)=V 2 q(l), W lq (l) = W 2 q{l) , 

w' lq (l) - Ul q (l) = P[w' 2q (l) - U2q{l)\, 

! n\ Vl i(0 _ a L m u 2 g (0 u 2g (0" 

t l -P\v 2 ,(<) z z 

We then using the general solution in Eq. (26), with coefficients A-|, .B 1; Ci for fluid #1, ensuring the correct decay 
in the far field, and coefficients A 2 , -B 2 , C 2 , D, E, F for fluid #2. 

Substituting this general solution into the boundary conditions in Eq. (52), we obtain a linear system with 9 
coefficients, 


(52a) 

(52b) 

(52c) 

(52d) 


pK 1 (p)A 2 + K 2 (p)B 2 + K 0 (p)C 2 + ph{p)D + I 2 (p)E + I 0 (p)F = 1, 
K 2 (p)B 2 - K 0 (p)C 2 + h{p)E - I 0 (p)F = - 1 , 
\pK 0 (p) - K\ (p)] A 2 + K 1 (p)B 2 + K!(p)C 2 - Iftlo(p) + h{p)] D 

-h{p)E - h{p)F = 0, 
lK\(l)Ai + K 2 (l)B 1 + KoiQC! = 0, 
IE\{V)A 2 + K 2 {l)B 2 + K 0 (l)C 2 + lI\{f)D + I 2 (l)E + /q(/)E = 0, 


K 2 (l)B 1 - i^o(0Ci - K 2 (l)B 2 + Kq[V)C 2 - O(0£ + MO* 1 = 0 , 
pATo(Z) - K \(/)] Ai + ^(/)Bi + Ad(0Ui - [Jtf o (0 - Ki(l)\A 2 - K x {l)B 2 - mi)C 2 

+ [lIo(l) + I\(l)]D + I\(l)E + I\(l)F = 0, 


2K 0 {1) - 2lK\(l) + 


~P 

+P 


2K 2 {1) - 21K\{1) + 
2Io(l) + 2lh(l)- 


Ki(D 

l 

mr 

1 

ui) 


A 1 


mo 

1 


-2K 2 (0 


a 2 — p 


mo 


- 2K 2 (0 


Bi- 

-P 


2K 0 (l) - 


mo 


2K o (0 


l 

mo 


1 


c\ 

C 2 


l 


D + p 


0(0 

l 


■ 2 / 2(0 


E + p 


2/0 (0 - 


0(0 

l 


F = 0 , 


-/fr(0^i - ^ 3 (0^i + Ad(0C7i + /?[Ad(0^2 + m0B 2 - K\(l)C 2 

+h(l)D - I 3 (l)E + h(l)F] = 0, 


(53a) 

(53b) 

(53c) 

(53d) 

(53e) 

(53f) 

(53g) 


(53h) 

(53i) 


which can be easily inverted numerically. 

The last step consists in moving to next order and computing the swimming speed. This is done similarly to the 
case with a finite slip length, and we apply Eq. (30) with A = 0. The average of w 2 ; on the filament is the swimming 
speed, U < ' 2 \ and we have 

u (2) = 1 

4 

which can be evaluated as 

C/ (2) {-B 2 K 0 (p) - C 2 K 0 (p) + A 2 [3K 0 (p) - pmp)\ 

—EIo(p) - FI 0 (p) - £>[3/ 0 (p) + ph(p)}} • 


W> 


2 q 


(p) + 


W2q(p) 


(54) 


(55) 
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(a) (b) 


FIG. 8: (Color online) Ratio between the swimming velocity in the two-fluid domain, U^, and that obtained in the case of a 
single Newtonian fluid, Uso', as a function of the mean dimensionless distance between the cylindrical filament and the fluid 
interface, kh\ (a): the dimensionless radius is fixed, kp = 0.1, and three values of viscosity ratios are considered (/3 = 0.01, 0.1 
and 0.5); (b): fixed viscosity ratio (/3 = 0.1), and three different values of dimensionless radius (kp = 0.01, 0.1 and 0.5). 


In Fig. 8 we plot the ratio between the swimming speed of the filament in the two-fluid domain, U^ 2 \ and that 
obtained in the case of a simple fluid, for a range of values of the dimensionless radius (kp, i.e. p) and viscosity 
ratio ((3). The results are reminiscent of the ones shown in Fig. 6 for the two-dimensional case. The swimming speed 
is always increased by the presence of a second fluid, potentially by order of magnitude when the wavelength of the 
swimmer and its radius are large compared to the thickness of the low-viscosity layer. 


IV. DISCUSSION 

In this paper, we presented a physical mechanism for the locomotion enhancement of microscopic swimmers in 
a complex fluid. The physical idea is that phase-separation near the surface of the swimmer leads to the presence 
of a low viscosity layer which promotes slip and decreases viscous friction. As a way to intuitively rationalize the 
results in our paper, we note - as is well known - that the locomotion in a fluid is governed by the ratio between the 
drag coefficients for the motion relative to the fluid perpendicular to and along the surface of the swimmer [6]. The 
presence of a fluidic interface, or of a finite slip length, affects normal hydrodynamic forces only weakly but leads to a 
systematic decrease of tangential viscous forces, and hence should always lead to faster swimming, as observed here. 
In a different context, but with some physical similarities, swimmers also always enhance their swimming speed in a 
network of polymer molecules [38] and in a porous medium [39]. 

Beyond the traditional geometrical assumptions made in our paper which are similar to a number of classical studies 
(namely solving the swimming problems for small amplitude motion and perfectly sinusoidal waveforms), one severe 
restriction of our two-fluid approach is the assumption that the interface between the two fluids remains flat. This is 
akin to saying that the time scale of the waving motion is much faster than the time scale for the readjustment of the 
interface, which is a reasonable assumption only for large fluid viscosities. That flat interface then provides an effective 
confinement to the swimmer, which is known to enhance locomotion [40]. A more sophisticated physicochemical model 
including molecular details of the phase separation would be required to solve for the dynamics of the thin film and 
to untangle the relative importance of viscosity difference vs. confinement in the increase of the swimming speed. 

What are the quantitative predictions of our models? In two dimensions, we obtained that the speed of a two- 
dimensional infinite swimming sheet is increased by 1 + 2/cA for the swimming with wave number k. Slip lengths 
of polymer solutions, A, have been measured in the range 0.1 — 10 pm [41]. For a microswimmer with wavelength 
A = 2n/k « 10 pm [6], the swimming speed can then be increased by 0(10%) up to by one order of magnitude. 
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In the three-dimensional case, we also obtained the increased speed shown in Fig. 4, which is consistent with the 
two-dimensional situation. Considering the same wavelength, and for a filament with radius p = 150 mn (so with 
dimensionless radius kp « 0.1), the predicted enhancement ranges from 0(30%) up to by forty times of the speed in 
the Newtonian fluid. 

In our second model, we used a two-fluids domain to describe wall depletion. In both two- and three-dimensions, 
we saw that when the wall depletion layer is very thin, the enhancement can be very large. A recent experiment 
by Gagnon, Shen and Arratia [42] considered the locomotion of the nematode C. elegans in concentrated polymer 
solutions and showed that the swimming speed can be increased significantly, by up to 65%. Is our model consistent 
with this result? The thickness of the low-viscosity layer in concentrated polymer solutions is a complex function of 
mean diameter of the particles and the concentration but can be estimated using the empirical formula [43] 



(56) 


where D p is the mean diameter of the particles, (j) the particle volume fraction, and 4> m the maximum packing fraction. 
The diameter of the particle (Xanthan gum in the experiment of Ref. [42]) is about 200 pm. For semi-concentrated 
solution, (j>/<j>m is 0.4, while for a concentrated solution the value is 0.8, so the slip layer thickness is around 120 /rrn 
and 40 pm respectively. The wavelength of the swimming worm is 2tt/ k ss 1 mm, and the dimensionless diameter kp 
is about 0.1. When the ratio of the viscosities /3 changes from 0.01 to 0.5, the enhancement predicted by our model 
is about 90 — 98%, which is less than a factor of two away from the experimental results, and indicates that our 
simplified approach captures the essential physics of the swimming enhancement. 
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